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Abstract. I review the current status of hadronic structure computations on the lattice. I describe 
the basic lattice techniques and difficulties and present some of the latest lattice results; in particular 
recent results of the RBC group using domain wall fermions are also discussed. 

Understanding the basie properties of matter requires the understanding of the nueleon 
strueture. Quantum Chromodynamies (QCD) is the theory deseribing strong interaetions 
and henee is responsible for the properties of the nuclear matter. Although QCD has 
been around for more than twenty years, its non-perturbative nature is an obstacle to the 
direct connection of low energy physics to quarks and gluons, the fundamental degrees 
of freedom of the theory. Unlike QED, non-perturbative techniques had to be developed 
in order to understand the QCD predictions at low energies. The lattice formulation of 
QCD is both a non-perturbative way to define the theory and a very powerful tool in 
understanding its properties. 

Deep inelastic scattering of leptons on nucleons has been an important tool in un¬ 
derstanding the structure of hadrons. Over the last few decades experiments at SLAC, 
Fermilab, CERN, DESY, and more recently at RHIC and JEAB, have measured the 
quark and gluon light cone distribution functions of the nucleon. These experiments 
have substantially advanced our knowledge of the properties of hadrons. However, we 
would also like to study how this observed rich phenomenology arises form first prin¬ 
ciples, i.e. QCD. With todays advances in computer technology, algorithms, and recent 
developments in lattice regularization of fermions, lattice calculations can complement 
the experimental effort and promote our understanding of the non-perturbative nature of 
QCD. 

The Lattice Formulation 

The continuum Euclidean path integral can be defined using the lattice regulator [HI] . 
In order to preserve gauge invariance the lattice gauge fields are link variables 

U^(;c) (1) 

where are the continuum gauge fields. The fermion fields live on the sites of the 
lattice. Naive discretization of the continuum fermionic action leads to the so-called 
fermion doubling problem. This problem can be avoided by either reinterpreting the 
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FIGURE 1. Deep inelastic scattering. 


additional light fermions as extra flavors (the Kogut-Susskind approaeh) or by intro¬ 
ducing an irrelevant dimension 5 operator that breaks chiral symmetry on the lattice 
and gives mass proportional to the inverse lattice cutoff to the fermion doublers (the 
Wilson approach). Recently, new lattice fermionic actions that both preserve chiral sym¬ 
metry on the lattice and do not suffer from the fermion doubling problem have been 
introduced. Such fermionic actions are the domain wall fermions the over¬ 

lap fermions [|^], and the fixed point fermions ^ . Having defined the lattice theory, 
correlation functions can be evaluated using Monte-Carlo integration in Euclidean space. 

However, parton distribution functions are defined in the Minkowski space, and hence 
cannot be directly computed in lattice QCD. Using the operator product expansion we 
can relate moments of the structure functions to forward matrix elements of gauge 
invariant local operators (for a pedagogical review see my These matrix elements 
can then be computed using lattice QCD. 

In a deep inelastic process (see Fig. [T|) the cross section is given by 


(f-o 

dQ.dE' 


1 a^E' 
q'^ E 




( 2 ) 


where is the lepton tensor, W^y is the hadronic tensor, q is the momentum transfer, 
is the nucleon mass. The initial and final energy and momentum of the lepton are 
{E,k) and {E',k') respectively. 

The hadronic tensor can be decomposed in the symmetric and anti-symmetric 

pieces: 




(3) 


The symmetric piece defines the unpolarized structure functions E^ and E 2 {E^ also 
for neutrino scattering). 


(x, Q^) - + ^^E,ix,Q^) + 
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while the anti-symmetric defines the polarized structure functions gj and g 2 


(4) 


=i£^''P'^qp (^(gi(x, Q^) +g2(x, Q^)) - ^ '^ g2(^. 6^)) ■ (5) 









where Pf^ and 5 ^ are the nucleon momentum and spin vectors, v = q-p, 

;c = 2^/2V and ~ —q^. 

At the leading twist in the operator product expansion the moments of the structure 
functions can be factorized at a scale /i in hard perturbative contributions (the Wilson 
coefficients) and low energy matrix elements of local gauge invariant operators: 



q=u,d 

q=u,d 

q=u,d 

' q=u,d 


where are the Wilson coefficients and are the non- 

perturbative matrix elements. At the leading twist and are related to the 

parton model distribution functions {x'^)q and 

(^-'>, = vW (7) 

In order to extract and we need to compute non-perturbatively 

the following matrix elements: 

i(*»pp,-*P,p»)pp2'"Pt..+'"-»xi 

( 8 ) 

{} implies symmetrization and [] anti-symmetrization of indices. The nucleon states 
Ip, 5 ) are normalized so that {p,s\p',s') = (2;r)^2£’(p)5(p — p')5^ ^ and — —m^. The 
operators are 


fi \^ ^ 

dY^,D^^---D^„q-trace 

f i\" _ ^ ^ 

I 2 j dYoYs ^M2 ■ ■ ■ d - trace 




( 9 ) 



where D= D — D and D, D are eovariant derivatives aeting on the right and the left 
respeetively. 

In Drell-Yan proeesses the transversity distribution ean be measured (for details 
see [ITTi [T3|] ). The relevant matrix element is 


and the operators 


are 





qy^Opy Dp^ - Dp„ q - trace. 


( 11 ) 


Lattice matrix elements 

In order to ealeulate on the lattiee the needed matrix elements we have to eompute 
nueleon three point funetions 

<x,P 


and nueleon two point funetions 


Ciptip^^) = E 

(x.p 


{Jn{p,t)Ja{p,0)) 

a,P 


(13) 


where J(p,0) and J(p,t) are ereation and annihilation operators of states with the 
quantum numbers of the nueleon. For unpolarized matrix elements F = while for 

the polarized F = —^* 75 7^ (k 7 ^ 4). The faetor is for projeeting out the positive 
parity part of the baryon propagator i.e. the nueleon. For the proton a typieal ehoiee is 


Ja(P,f) 


E » 


ip-x^abc 


u^(x,t)Cy^d^(x,t) 


u%(x,t) 


(14) 


where C = 7472 the eharge eonjugation matrix, a is a spinor index and a,b,c are eolor 
indiees. When tT 3 > 0 


c ( 3 ,) - Z + 

'-'ipppq) ~ E{p) ^ ^ 

C[pf(p,FT) = Y,'^ap^NiP^s){P^s\^\P^s)U^ip.s)e~^Nipy + ■■■ (15) 

OC,p,S 


where U{p,s) is the nueleon spinor whieh satisfies the Dirae equation and 
{Q\Ja{pd)\p^s) = y/Z^U‘^{p,s). From Eq. [1^ and Eq. |] (or Eq. the required 
matrix elements ean be extraeted from the ratio of three point funetions over two point 










functions. In practice we would like to aehieve the asymptotie behavior of Eq. with as 
small as possible t and x. For that reason the interpolating operator J is tuned so that the 
overlap with the exited nucleon states would be as small as possible. For more details 
on the teehnieal aspeets of the lattiee ealeulation the reader may refer to [|T^, [T^ . 

In order to reduee the eomputational eost of ealeulating the above eorrelation fune- 
tions some times the so-ealled quenched approximation is used. In this approximation 
the quark loop contributions to the path integral are ignored. Quenehing reduees the 
eomputational eost by several orders of magnitude, while for eertain quantities it in- 
troduees a systematie error ~ 10%. [| In addition, lattiee eomputations are typieally per¬ 
formed with heavier quark masses than the physieal up and down quarks. Henee we have 
to perform extrapolations to the chiral limit. If the quark masses are light enough, ehiral 
perturbation theory [ ]T8| , [T^ , ^ ean be used to calculate the dependence of the matrix el¬ 
ements on the quark massj^ Finally the lattiee matrix elements have to be renormalized, 
typieally to MS, and extrapolated to the eontinuum limit. 


Renormalization 


The renormalized operators at scale n are obtained from the lattice operators ealeu- 
lated at lattiee spaeing a from 


(16) 

in the ease of multiplieatively renormalized operators. In general, there is operator 
mixing and as a result the above relation becomes 


Xm) =z(m;«) 


slat 
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(17) 


where are a set of operators allowed by symmetries to mix, and is the dimension 
of eaeh operator. It is evident that if mixing with lower dimensional operators occur, the 
mixing eoeffieients are power divergent as we approaeh the continuum limit. Hence we 
have to eompute these terms non-perturbatively in order to aeeurately renormalize the 
operators. Higher dimensional operators are typieally ignored sinee their effeets vanish 
in the continuum limit. In certain eases we may want to compute these eoeffieients in 
order to remove part of the systematic error introduced by the finite eutoff. 

The mixing of lattiee operators is more eomplieated than that of the eontinuum 
operators, sinee on the lattiee we do not have all the eontinuum symmetries. In partieular, 
0(4) rotational symmetry in Euclidean spaee is broken down to the hypereubie group 
H{A). As a result, an irredueible representation of 0(4) is redueible under H{A) and 
henee mixing of operators that would not oeeur in the eontinuum ean oeeur on the lattiee. 


For a detailed analysis of the 7/(4) group representations see [0, ^ and referenees 
therein. In lattice ealeulations we have to seleet carefully the lattice operators so that 


^ Note that there are quantities for which the quenched approximation introduces uncontrollable errors. 
^ In the case of the quenched approximation the so-called quenched chiral permrbation theory is used. 
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FIGURE 2 . Quark density {x)q vs. the pion mass squared, [left] The connected up (octagons) and down 
(diamonds) quark contributions, [right] The flavor non-singlet 


mixing with lower dimensional operators does not occur and hence no power divergent 
coefficients in Eq. |T^ are encountered. This turns out to be a significant constraint on 
how many moments can be practically computed on the lattice. 

Another symmetry that is broken on the lattice for Wilson fermions is chiral symme¬ 
try. This results in mixings with lower dimensional operators for the dn matrix elements. 
Fortunately, in this case we can use lattice fermions, such as domain wall or overlap and 
fixed point fermions, that respect chiral symmetry on the lattice. For Wilson fermions, 
the renormalization of ^2 has been done non-perturbatively as described in [ Jlby . 

The renormalization constants for all the operators relevant to structure function cal¬ 
culations have been computed perturbatively for Wilson fermions, improved and unim¬ 
proved |2^]. Moreover, the RI-MOM scheme has been used to renormalize non- 

perturbatively both local and derivative operators [ ^7] , In the Schroedinger 
functional scheme (developed by the AFPHA collaboration), all local operator renor¬ 
malizations and the renormalization of V 2 have been computed 0B- In addition, 
work is underway for computing the constants for flavor singlet operators [|^]. For do¬ 
main wall fermions, all local operators have been renormalized non-perturbatively [ ]77| ] 
using the RI-MOM scheme, and also perturbatively 


LATTICE RESULTS 


In the last several years, the lattice community (QCDSF/UKQCD and FHP/SESAM col¬ 
laborations) has made a substantial effort to compute the first few moments of the nu¬ 
cleon structure functions. Apart from the constraints imposed by the renormalization of 
the operators mentioned above, the requirement of having nucleon states with non-zero 
momentum^ limits the number of moments we can compute. These are the first three 
moments of the unpolarized structure functions, the first two moments of the polarized 
structure functions, and the first two moments of the transversity. These computations 
have been performed both in quenched and in full QCD with improved and unimproved 
Wilson fermions [|T^, [1^, . 

The RBC group has recently begun quenched computations with domain wall 


fermions [p5|] . Our current results are restricted only to those matrix elements that can 


^ operators with more than one derivative need non-zero momentum nucleon states see Eq. ^ and Eq. 











be computed with zero momentum nucleon states. We use the DBW2 gauge action 
which is known to improve the domain wall fermion chiral properties ^ . We have 
416 lattices of size 16^ x 32 at /3 = 0.870 with lattice spacing = 1.3GeV, providing 
us with a physical volume (~ (2.4/m)large enough to reduce finite size effects known 
to affect some nucleon matrix elements, such as Using fifth dimension 

length = 16 we achieve a residual mass mres ~ O.SMeV [|3^, The input quark 
masses ranged from 0.02 to 0.10, providing pion masses ranging from 390MeV to 


850MeV. Further technical details of our calculation can be found in [35]. 


Unpolarized Structure Functions 


The first three moments of the unpolarized structure functions have been computed 
by QCDSF in the quenched approximation. The needed chiral and continuum extrapo¬ 
lations have also been performed. A summary of recent results can be found in pOtj. In 
comparison with MRS phenomenological results, the lattice results are typically higher. 
Also, V 3 is smaller than V 4 , while V 3 > V 4 is phenomenologically expected. The same 
computations have been performed by LHP/SES AM in full QCD [ p^ . Their results in¬ 
dicate that dynamical fermions have only a small effect on the matrix elements they 
studied. 

It has been argued that the main reason for such discrepancies is the fact that lattice 
computations are performed at rather heavy quark masses and then extrapolated linearly 
to the chiral limit [^H, @ 1 , p3|] . For that reason, we need computations at much lighter 
quark masses in order to see whether there is a disagreement with phenomenological 
expectations. In quenched QCD, a study with very light quark masses has been done [ ]44| ] 
indicating that the linear behavior persists down to 300MeV pion masses. 

In Fig. ^ we present our results for the quark density distribution {x)q (V2). We plot the 
unrenormalized result for {x)u, {x)^ and the flavor non-singlet {x)^_^. Down to 380MeV 
pion mass no significant curvature within our statistical errors can be seen.[| The ratio 
{x)u/ {x)^ is 2.41(4), linearly extrapolated to the chiral limit, is in agreement with the 
quenched Wilson fermion results [ P3| , p^ . 

Polarized Structure Functions 


The nucleon axial charge is related to the first moment of the polarized structure 
function gj. The current experimental value for gj^^/gy measured from neutron beta 
decays is 1.2670(30) [^3]]. Lattice calculations, quenched and dynamical, have been 
underestimating this quantity typically by 10% to 20% [T^, ^ [T^, |T^ ^]. For 

earlier calculations see also 

One of the systematic errors believed to affect these calculations is the finite volume. 
In order to study this effect we performed two calculations. One with spatial volume 
2.4^/m^ and another with spatial volume 1.2^/m^. Our results are shown in Fig. ^[right] . 
Between these two volumes it is clear that there is a finite volume effect of about 20% at 


In [ 361 we had an indication of some curvature but this effect went away as we doubled the statistics. 
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FIGURE 3. Helicity (l)^^ vs. the pion mass squared, [left] The connected up (octagons) and down 
(diamonds) quark contributions, [right] The nucleon axial charge i.e. flavor non-singlet 
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FIGURE 4. Helicity (x)^^ vs. the pion mass squared, [left] The connected up (octagons) and down 
(diamonds) quark contributions, [right] The flavor non-singlet 


the chiral limit. In addition, the linearly extrapolated to the chiral limit value for gp^/gy 
is 1.21(2). For a detailed analysis of this computation see [^. Note that for domain 
wall fermions gjgy does not require renormalization, since the finite renormalization 
constants of the axial and the vector currents Z^, Zy are equal []32|, In Fig. ||[left] we 
present the up and down quark contributions of (1)^^ for the proton renormalized using 
= .77759(45) [lJ7l] . In Fig. ^ we present our unrenormalized data for The ratio 
m linearly extrapolated to the chiral limit is roughly —4, consistent with other 
lattice results MM- The lowest moment of the transversity (1)^^ is also measured. In 
Fig. 1^ we plot the unrenormalized contributions for both the up and down quark, and the 
flavor non-singlet combination Again the quark mass dependence is very mild 

and there is no sign of a chiral log behavior. The ratio linearly extrapolated 

to the chiral limit is also roughly —4. 

For computing moments of ^2 we need to calculate the twist 3 matrix elements dn- 
We computed the matrix element which contributes to the first moment of ^ 2 - If 


chiral symmetry is broken the operator 


34 
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q which is used to 


compute d^ mixes with the lower dimensional operator Hence in Wilson 

fermion calculations a non perturbative subtraction has to be performed. This has been 
done for ^2 by QCDSF [[T^, With domain wall fermions this kind of mixing is 
proportional to the residual mass (~ mres/^i), which in our case is negligible. Thus we 
expect that a straightforward computation of d^ with domain wall fermions provides 
directly the physically interesting result. In Fig. § we present our unrenormalized results 
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FIGURE 5. Transversity (l)^^ vs. the pion mass squared, [left] The connected up (octagons) and down 
(diamonds) quark contributions, [right] The flavor non-singlet 
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for as a function of the quark mass. For comparison we also plot the unsubtracted 
quenched Wilson results for /3 = 6.0 from [|T7|]. The fact that our result almost vanishes 
at the chiral limit is an indication that the power divergent mixing is absent for domain 
wall fermions. The behavior we find for the matrix element is consistent with that of 
the subtracted ^2 computed by QCDSF [jl^ ^ with Wilson fermions. 

CONCLUSIONS 


In conclusion, lattice computations can play an important role in understanding the 
hadronic structure and the fundamental properties of QCD. Although some difficulties 
still exist, several significant steps have been made. Advances in computer technology 
are expected to play a significant role in pushing these computations closer to the chiral 
limit and in including dynamical fermions. RBC has already begun preliminary dynam¬ 
ical domain wall fermion computations [ ]5T1 ] which we expect to be pushed forward with 
the arrival of QCDOC In the near future, we also expect to complete the non- 
perturbative renormalization of the relevant derivative operators in quenched QCD. 
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